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Baryon number cumulants are invaluable tools to diagnose the primordial stage of heavy ion 
collisions if they can be measured. In experiments, however, proton number cumulants have been 
measured as substitutes. In fact, proton number fluctuations are further modified in the hadron 
phase and are different from those of the baryon number. We show that the isospin distribution of 
nucleons at kinetic freeze-out is binomial and factorized. This leads to formulas that express the 
baryon number cumulants solely in terms of proton number fluctuations, which are experimentally 
observable. 
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The order of the phase transition of quantum chromo- 
dynamics (QCD) at nonzero temperature (T) is believed 
to change from crossover [ij to first order at a nonzero 
baryon chemical potential (hb)- The existence of the 
QCD critical point is thus expected in the phase dia- 
gram on the T-/iB plane [4I. Experiments to explore the 
phase structure at nonzero /^b, especially the existence 
of the critical point, are now ongoing in the energy scan 
program at the Relativistic Heavy Ion Collider (RHIC) 
[3, 3] , and will also be performed in future facilities [5|, Q . 
Much attention has also been paid to thisproblem from 
numerical experiments on the lattice [l|, lZ| . The estab- 
lishment of the QCD phase structure at nonzero /ib is 
an important issue, not only to deepen our knowledge of 
the matter described by QCD, but also to gain under- 
standing of a wide array of topics in physics which share 
the concepts of phase transitions and techniques to treat 
strongly correlated many-body systems. 

Fluctuations, which are experimentally measured by 
event-by-event analyses in heavy ion collisions, are 
promising observables to probe the properties of created 
fireballs [SJ], as their behaviors are sensitive to the state 
of the matter. For example, because of the singularity at 
the critical point, fluctuations of various physical quan- 
tities, including skewness and kurtosis, behave anoma- 
lously near the critical point [9l4ll|. One can also argue 
that ratios between the cumulants of conserved charges 
are sensitive to the magnitudes of the charge carried 
by the quasiparticles composing the system, and hence 
they behave differently in the hadronic and quark-gluon 
phases [134l4|. Recently, it was also pointed out that 
some higher-order cumulants of conserved charges change 
signs around the phase boundary of QCD, which would 
serve as clear experimental signatures to determine the 
location of the matter in the phase diagram 15l4l7|. 

Among the fluctuation observables, those of conserved 
charges can reflect fluctuations produced in earlier stages 
during the time evolution of fireballs, than non-conserved 
ones [l8|. This is because the variation of a conserved 
charge in a volume is achieved only through diffusion. 



which makes the relaxation to equilibrium slower. In 
fact, it is argued that if the rapidity range of a detector 
is taken to be sufficiently large, whereas the range should 
be kept narrow enough so that the rest of fireballs can be 
regarded as the heat bath, the effects of diffusion are well 
suppressed and fluctuations produced in t he q uark-gluon 
phase can be detected experimentally |12l.ll3|. 

The dependences of the proton number fluctuations, 
cumulants up to fourth order, on the beam energy ^/s, 
have been recently measured by the STAR collabora- 
tion at RHIC [3|, y. The result appears to be almost 
consistent with the prediction of the hadron resonance 
gas (HRG) model jl9[; although the experimental result 
shows some deviation from the prediction at small ^/s, 
it is at most of the order of 20% [J|. The proton num- 
ber, however, is not a conserved quantity, and in fact 
we will see later that its fluctuations significantly evolve 
in the hadronic stage, which makes the experimentally 
measured fluctuations close to those in the equilibrated 
hadronic matter. The agreement between the experi- 
ments and the HRG model in the proton number fluctu- 
ations I4J is in part due to these effects, and hence it does 
not immediately exclude the slow baryon number diffu- 
sion in the hadronic stage. Although the measurement of 
the baryon number, which is a conserved charge, is desir- 
able to probe fluctuations generated in earlier stages, its 
direct experimental measurement has been considered to 
be impossible because of the difficulty in detecting and 
identifying neutrons. 

In this Rapid Communication, we show that the exper- 
imentally measured proton number fluctuations are nev- 
ertheless directly related to baryon number fluctuations 
in earlier stages, and we present concrete formulas that 
relate the baryon number cumulants and these experi- 
mental observables. The key observation is that the dis- 
tributions of (anti-)proton and (anti-)neutron numbers 
in the final state are well described by binomial distribu- 
tions. As will be argued in detail later, this observation 
is well justified at least for RHIC energy, and is expected 
to hold for Vi > lOGeV. 



Experimentally, the electric charge can be measured 
directly. Electric charge fluctuations, however, contain 
the contribution of isospin fluctuations, which are non- 
singular at the critical point, in addition to baryon num- 
ber fluctuations [10| . The signals of the phase transition 
in this observable thus generally become weak owing to 
the non-singular contribution (such a tendency, for ex- 
ample, in the third moments is seen in Ref. [la]). In this 
sense, the baryon number fluctuations are superior to the 
electric ones as probes of the QCD phase structure. 

Throughout this Rapid Communication, we use Nx to 
represent the number of particles X leaving the system 
after each collision event, where X = p, n, and B repre- 
sent proton, neutron, and baryon, respectively, and their 
anti-particles, p, n, and B. The net and total numbers 
are denoted as n'^x"^'^ = Nx^N^ and iv]f*°*^ = Nx+Nx, 
respectively. 

Before starting the main discussion on the cumulants of 
baryon and proton numbers, let us briefly consider how 
the proton number fluctuations evolve in the hadronic 
stage. The most important process responsible for the 
variation of the proton number is the charge exchange 
reactions with thermal pious mediated by A"*" (1232) and 
A°(1232) resonances: 
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Because of the small energy required and the large cross 
sections, these reactions proceed even after chemical 
freeze-out, as we demonstrate later. We note that these 
reactions do not alter the average abundances {Np) and 
{Np) if the isospin chemical potential vanishes, while 
they modify the fluctuations of Np and Np. Because 
chemical freeze-out is a concept that describes ratios be- 
tween particle abundances such as {Np)/{Np), these re- 
actions below the chemical freeze-out temperature Tchcm 
do not contradict the statistical model. The success of 
the model, on the other hand, indicates that the creation 
and annihilation of (anti-)nucleons hardly occur below 

-^ chem- 

The importance of reactions ([!]) below Tchcm is con- 
flrmed by evaluating the mean time of the nucleons for 
these reactions. Provided that the pious have a thermal 
distribution, the mean time r that a proton at rest in the 
medium forms A"*" or A", being scattered by a thermal 
pion, is evaluated to be 



(27r)3 
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with the Bose distribution function n{E) — {e^^"^ — 1)^^, 
pion velocity v-^ = fc^/^E^, E^, ~ a/to^ -|- k1, and the pion 
mass TO^. a{Ec.n\) is the sum of the cross sections for 
pi: reactions producing A+ and A*^ with a center-of-mass 
energy -Ec.m. = [("iN + E^^Y ^ ^'iV'^'^ ^^lih the nucleon 
mass TTiN. To evaluate Eq. ^, we assume a cross section 
of Breit-Wigner type, cj{E,,^) = CTA(rV4)/((^c.m. - 



E/^Y +r^/4), which is a sufflcient approximation for our 
purpose, with hadron properties in the vacuum, m^ — 
940MeV, m^ = 140MeV, E/^ = 1232McV, T = llOMeV, 
and cta = 20fm^ [20|. The mean time is then evaluated 
to be 3 - 4 fm for T = 150 - 170MeV. One can also 
check that this mean time hardly changes even for moving 
protons in the range of momentum p < 3T. On the 
other hand, dynamical models for RHIC energy predict 
that protons stay in the hadronic gas and continue to 
interact for several tens of fm on average at midrapidity 



21| . which is significantly longer than the mean time 
and the lifetime of A, l/F ~ 1.8fm. This result shows 
that nucleons in the fireball indeed undergo this reaction 
several times on average in the hadronic stage [2JJ. The 
ratio of the probabilities that a proton in the medium 
produces a A+ or A*' and then decays into p and n is 
5 : 4, which is determined by the isospin SU(2) algebra. 
Whereas this probability is not even, after repeating the 
above processes several times in the hadronic stage, the 
nucleons tend to completely forget their initial isospin. 

The above discussion shows that the evolution of pro- 
ton number fiuctuations in the hadronic stage is domi- 
nantly made via the exchanges of the two isospin states 
of the nucleons. Now, we further assert that isospins 
of all nucleons in the final state are uncorrelated. This 
statement is well justified when the hadronic medium 
fulfills the following two conditions: (i) The medium ef- 
fects on the branching ratios and formation rates of A 
are insensitive to the proton and neutron number densi- 
ties Up and n„ (and the same holds for the anti-particle 
sector as well), and (ii) (anti-)nucleon-(anti-)nucleon in- 
teractions generating correlations between two nucleons 
hardly occur. As we will see later, these two conditions 
are well satisfied below Tchcm except for low-energy colli- 
sions. The probability distribution of finding Np and Nn 
{Np and Nn) particles in the final state in each event 
then becomes binomial. Under the isospin symmetry 
[19j, this fact enables to factorize the probability distri- 
bution P{Np,Nn,Np,Nn) having Np, iV„, Np, and Nn 
particles in each event as 



P{Np,Nn,Np,N^) 
- F{NB,Nj,)B{Np; Nb)B{Np; iVg), 



where B{k;N) = 2-^N\/{kl{N - k)\) is the 
distribution function with an equal probability, 
right-hand side (RHS) of Eq. ([3]) we have used 
Nq defined by A^^b = Np+Nn and N^ = Np+Nf,. 
later elucidate this notation to use the baryon 
Nb and A^g in place of the nucleon numbers 
the probability distribution Eq. (|3]), the event- 
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average of a function f{Np, Np) is given by 

{f{Np,Np))^ Y, PiNp,N„,Np,N^)fiNp,Np) 

^{p.n.P.K} 

= J2 F{Nb,N^) J2 f(Np.Np) 

xBiNp;NB)BiNp;N^). (4) 

The factorization Eq. ([3]) leads to 

{Nt'^)-l{Nt'^), (5) 

^(^^(„et))2^ J_l^^SNt'^f) + \{NJ^°'), (6) 

{{SNJ,'-^^f) =i((5<-))3) + ^(5Ar^-)<5iV^-)), (7) 
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+ i^{mt'^r)-UNt'^), (8) 



where (SiVx ^ Nx ~ {Nx}. To derive Eqs. ©-([5]), we 
have used the fact that the sums over Np and Np in 
Eq. (j31) can be taken separately with corresponding bi- 
nomial functions, e.g., J^n ^pB{Np;Nb) = Nb/2 and 
Eat^ N^BiNp; Nb) = iV|/4 + Nb/4. 

Equation ([3]) also enables to represent the baryon num- 
ber cumulants by those of the net and total proton num- 
bers as 
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where we have used relations for mixed cumulants such 
as {SN^'''''>dN^°''>) = 4{SNJr"'''SN^'°'^)-2{N^'°'^) which 
are obtained with Eq. ^. Since the RHSs of Eqs. dH) 
- (fT2|) consist of only Np and Np ° , which are ex- 
perimentally observable, these are formulas that express 
baryon number cumulants solely in terms of experimen- 
tal observables. We remind that no specific form of 
F{Nb, Nb) is assumed in deriving these results. 

We remark that N^''^^ (iV^*°*^) in Eqs. © - ([121) are 
interpreted to be the sum of all net (total) baryon num- 
bers entering a region in the phase space in the final state 
of each event. If the diffusion of the baryon number in 
the hadronic stage is slow [l^, HSj , the information on the 



primordial fluctuations remains in F{Nb, Nb) in Eq. ^ 
and, as a result, in baryon number cumulants. 

Next, let us inspect the validity of Eq. ([3]) in more 
detail. First, we consider the conditions (i) and (ii) in- 
troduced above Eq. ([3]). In the medium, the decay rate 
of A acquires the statistical factor 



il-f{EN)){l + n{E^)), 



(13) 



where f{E) = {e^E-t,B)/T _^ ^yi jg ^-^^ Fermi distribu- 
tion function and E^ and E^^ are the energies of the 
nucleon and pion produced by the decay, respectively. 
The first term in Eq. P^ represents the Pauli blocking 
effect. At RHIC energy, the Boltzmann approximation 
is well applied to nucleons below Tchcm since T <^ mpf 
and |/^b| ^ m-N- Thus, the Pauli blocking effect can be 
almost ignored. The Bose factor (1 -|-n(£',r)) in Eq. (fT3l) . 
on the other hand, has a nonnegligible contribution since 
TOtt — Tchom- The density of the pions, however, is more 
than one order larger than that of the nucleons below 
Tchcm- The Bose factor thus must be insensitive to Up 
and n„, while it leads to the enhancement of the decay 
of A in the medium, which acts in favor of the isospin 
randomization. The large pion density also means that 
the mean time for a nucleon to form A is insensitive to Up 
and n„. Condition (i) is thus well satisfied below Tchom at 
RHIC energy. The validity of condition (ii) is conjectured 
from the success of the statistical model as follows. The 
statistical model indicates that the pair annihilation of a 
N and an N terminates at Tchcm- NN and NN reactions 
are then also expected to terminate there, because the 
elastic cross section of NN is significantly smaller than 
the inelastic one, and the total cross section of NN be- 
haves similarly to that of NN for ^c.m. < IGeV ^. 
Condition (ii) thus should also be satisfied for T < Tchcm - 
Intuitively speaking, in a hot medium the nucleons are so 
dilutely distributed that they do not feel one another's 
existence, while there are so many pions which can be 
regarded as the heat bath when the nucleon sector is 
concerned. The large pion density also enables to use 
the binomial distribution independently of the initial nu- 
cleon isospin density. 

Second, while so far we have limited our attention to 
the nucleon reactions mediated by A, other interactions 
can also take place in the medium. It is also possible 
that A interacts with a thermal pion to form another 
resonance before the decay [24I . All these reactions with 
thermal pions, however, proceed with a certain probabil- 
ity determined by the isospin SU(2) symmetry as long 
as they are caused by the strong interaction, and the 
reactions of a baryon make its isospin random. Strange 
baryons, on the other hand, decay via the weak or electro- 
magnetic interaction outside the fireball. In particular, A 
and S are important among them. A decays into p and n 
with a branching ratio of 16 : 9. Provided that the three 
isospin states of E are produced with an equal probability 
in the medium, the ratio of probabilities that a E decays 



into p and n is about 1 : 1.6 [20|. Although these ratios 
are not even, because the abundances of A and S are 
smaU compared to the nucleons, to a first approximation 
it is suitable for our purpose to regard these probabili- 
ties to be equal and to incorporate nucleons produced by 
the decays of A and S in Np and iV„ in Eq. ([3]). This 
promotes the nucleon numbers to those of the baryons 
in Eq. ^. The treatment of strange baryons, however, 
may require more detailed arguments, especially on their 
quantitative effects on higher-order cumulants, which will 
be addressed elsewhere. Inclusion of higher baryonic res- 
onances and light nuclei such as deuterons will not affect 
our conclusions owing to their negligible abundances. 

While the factorization Eq. (j3|) is fully established for 
RHIC energy, the binomiality will eventually break down 
as the beam energy is decreased. At very low beam en- 
ergy, pions are not produced enough and nucleons will 
not undergo charge exchange reactions sufficiently below 
Tchcin- We deduce that this happens when Tchom ^ Jti^. 
When the reactions hardly occur, the isospin correlations 
generated at the hadronization will remain until the final 
state. At low beam energy, also the nucleon density be- 
comes comparable with that of the pions, and the latter 
can no longer be regarded as the heat bath to absorb the 
isospin fluctuations of the former. From the ^/s depen- 
dence of the chemical freeze-out line on the T-fi^ plane 
[23j, and considering the validity of these two conditions, 
we deduce that Eq. ^ is well applicable to the range of 
beam energy -^/s > 10 GeV. 

In the argument to derive Eqs. (|5|) — (fT2|) . we have im- 
plicitly assumed that the hadronic medium is isospin 
symmetric. While the effect of nonzero isospin density 
should be well suppressed for large -^i where a large 
number of particles having nonzero isospin charges are 
produced, at lower energies this effect gives rise to a non- 
negligible modification of Eqs. ([5|) — (fT2|) . When the sys- 
tem has nonzero isospin density, the probability that a 
nucleon at the early stage of the hadron phase becomes 
a proton or a neutron in the final state is no longer even. 
This effect is, as long as conditions (i) and (ii) introduced 
above Eq. ([3]) hold, incorporated into our results by sim- 
ply replacing the binomial function B{Np; Nb) in Eq. ([3]) 
with that having a probability k = {Np)/{Np + Nn), and 
a similar replacement to B{Np; N-g). Our explicit anal- 
ysis indicates that the effect of nonzero isospin density 
on Eqs. ([5]) — (fT2|) is relatively small and well suppressed 
when Tchcm > 'Ti^r and a sufficient number of pions hav- 
ing isospin charges are produced at chemical freeze-out. 
Since this modification requires a straightforward but 
lengthy calculation, we will elucidate the analysis in a 
forthcoming paper. 

Now, let us apply our results to the latest experimen- 
tal data from STAR [3, 4.]- To estimate how the bino- 
mial nature of nucleon isospins affects the proton number 
fluctuations, we first consider Eqs. (IS]) — ®. In order to 

.(tot) 



these equations, we temporarily postulate that Nb and 
iVg have thermal distributions fixed at chemical freeze- 
out as the statistical model suggests, while the distri- 
bution of their combination, Nj^'^ ', deviates from the 
thermal one reflecting the baryon number conservation. 
Under this assumption, the distributions of A'^b and Np 
are Poissonian, and hence the cumulants of the baryon 
and proton numbers satisfy 



{Nb) = {{SNb)') = {{SNb)') 

= 2{Np)iiG = 2{{SNpf)iiG = 2{{6Np)^) 



HG, 



(14) 



and the same for anti-baryon numbers, where (•)hg is the 
expectation value for free hadron gas (HG) composed of 
mesons and nucleons at Tchom, *-e- a simplified version 
of the HRG model jl3] . Equations (|6|) and ([7|) are then 
expressed as 



im^ 
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"°'^^')hg, (15) 



((5A.(not))3^ = ^((^A^^"¥) + ^((5A('-))3)hg. (16) 

To derive these results, we decomposed, for example, the 
second term in Eq. ^ as 

r(not)rAr(tot) 



- 2{{6Np)^)^G - 2{{SNpf)BG = 2((<5A^^-'))3) 



HG, 



(17) 

The results in Eqs. (fTS]) and ([T6| show that the second 
terms on the RHSs, which come from the binomial distri- 
butions of nucleon isospin, make a large contribution to 
the cumulants of the proton number, and they become 
more significant as the order increases. Although one 
cannot derive a similar result for the fourth-order rela- 
tion, from the factor 1/16 in the first term of Eq. ^ it 
is clear that the effect of the fourth-order baryon num- 
ber cumulant on the proton number one is more sup- 
pressed in this order. The suppression of the first term 
in Eqs. ([6]) — (H)) may be one of the reasons why the results 
of the STAR experiment and the HRG model appear to 
be consistent with each other. In this sense, it is interest- 
ing that the experimental results for skewness and kur- 
tosis have small but significant deviations from the HRG 
predictions at -^s < 50 GeV [J| . The deviation, for exam- 
ple, in skewness, can be a consequence of ((JA^g )'^) in 



Eq. ()16|) . which possibly refiects the properties of mat- 
ter in the early stage. Baryon number cumulants are, 
of course, directly determined with experimental observ- 
ables using Eqs. ([9]) — (1121 . It is worth emphasizing that 
the RHSs of Eqs. (fTTj) and (IT2|) have terms which would 
lead to the negative cumulants discussed in Refs. [l5l[l7|. 

or the suppression of the ratio ((5a4"''*V)/((<5^b"''*^)^) 



We note that when the distribution of A^, 



(not) 



also fol- 



estimate the contributions of terms including N^ in lows that in the HG, in addition to the above postulation. 



the RHSs of Eqs. ^ and dH]) reduce to ((5a4"''*^)^)hg 
and {{SNp )'^)hg- A way to check this is to use the 
fact that the (anti-)nucleon numbers in the HG are well 
described by the Poisson distribution owing to the Boltz- 
niann approximation, and that the Poisson distribution 
with an average A, P\{N), satisfies P\{Ni)P\{N2) — 
P2\{Ni + N2)B{Ni;Ni + N2). The HG thus corresponds 
to a special case of Eq. ([3]), where 



F(7VB,iVg) = P^.iNB)PiN^)iN^) 



(18) 



In this Rapid Communication, we derived relations 
between the baryon and proton number cumulants, 
Eqs. ©-(El) and ©-([ni), respectively, on the basis of 
the binomial nature of (anti-)nucleon isospin numbers in 
the final state. These results enable to immediately de- 
termine the baryon number cumulants with experimental 
results in heavy ion collisions, which will provide signifi- 
cant information about the QCD phase diagram. Though 
these results are obtained for the isosymmetric case, in- 
corporation of nonzero isospin density is straightforward 
and will be discussed elsewhere. 
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